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X-ray observations of the galaxy cluster Abell 2029 to the virial 
radius 
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ABSTRACT 

We present Suzaku observations of the galaxy cluster Abell 2029, which exploit Suzaku's low 
particle background to probe the ICM at radii beyond that which is possible with previous 
observations, and with better azimuthal coverage. We find significant anisotropics in the tem- 
perature and entropy profiles, with a region of lower temperature and entropy occurring to the 
south east, possibly the result of accretion activity in this direction. Away from this cold fea- 
ture, the thermodynamic properties are consistent with an entropy profile which rises, but less 
steeply than the predictions of purely gravitational hierarchical structure formation. Excess 
emission in the northern direction can be explained due to the overlap of the emission from 
the outskirts of Abell 2029 and nearby Abell 2033 (which is at slightly higher redshift). These 
observations suggest that the assumptions of spherical symmetry and hydrostatic equilibrium 
break down in the outskirts of galaxy clusters, which poses challenges for modelling clus- 
ter masses at large radii and presents opportunities for studying the formation and accretion 
history of clusters. 
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1 INTRODUCTION 

, Suzaku's ability to probe the intracluster medium (ICM) to the 

■ virial radius has been shown for a number of relaxed clusters 

■ so far taeorge et al."2009'; 'Reiprich et al."2009^, 'Bautz et al."2009'; 

jHoshino et al. 2010; Kawaharada et al. 2010; Simionescu et al. 

l2011 ii). The study of the outskirts of galaxy clusters has been pre- 
vented until recently due to its low surface brightness, and due to 
the high levels of background radiation to which X-ray telescopes 
are exposed. The virial radius broadly represents the boundary be- 
tween the central cluster which is in hydrostatic equilibrium, and 
the outskirts where matter is still accreting onto the cluster as it 
continues to form. These measurements are therefore crucial in un- 
derstanding the formation history of galaxy clusters, and in con- 
straining models which describe the formation process. They also 
allow mass estimates of galaxy clusters to be improved by allow- 
ing the dynamic state of the ICM to be understood to greater radii, 
whereas previous mass determinations have extrapolated measure- 
ments of the centres of clusters. This increases the accuracy with 
which galaxy cluster s can be used as prob es of cosmological pa - 
rameters (such as in lAUen et al.ll2008L and IVikhUnin et~Sll2009h . 
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Here we add Abell 2029 to the growing collection of clusters which 
have been investigated to the virial radius. 



Abell 2029 (z=0.0767) is a relaxed galaxy cluster which 
has been studied ex tensively in X-rays with ASCA and 
JSarazin et al.l JT998, hereafter 
ri999.). Chandra 



ROSAT JSarazin et al.l JTp' 
(iMolendi & De Grandil flgS 



Clarke et al .112004 



S98), 
fLewis 



Bep poSAX 
et al.l I2OO2I ; 



Vikhlininetal.11200 5. 2006), and XMM-Newton 
( ISnow denetal."2008; Bourdin & Mazzotta 2008*). It has a large 
cD galaxy ( Uson et al. 1991) whose major axis is aligned in the 
NE to SW direction, in approximately the same direction as that 
joining it to nearby Abell 2033. 



Along with Abell 2033, Abell 2028 and Abell 2066, 
Abell 2029 forms a small supercluster (catalogue number 154, 
lEinasto et alj|200in . the closest cluster of which is Abell 2033 at 
35 arcmin to the North, as shown in Fig. [T] Abell 2029 provides an 
excellent opportunity to study how the outskirts of a galaxy clus- 
ter are affected by its surroundings. Gastaldello et al. (2010) has 
suggested that a further cluster, IRXS J150935.9+084605, is asso- 
ciated with the supercluster, and their locations are shown on Fig. 

m 
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Figure 1. Suzaku pointings of A2029 overlaid on a ROSAT PSPC image 
of Abell 2029. Tlie five pointings are in green, the 6tli pointing from tlie 
public arcliive is in red. Abell 2033 is highlighted by the blue circle to the 
North. The blue annuli shown are used in extracting spectra and have ring 
radii at 2.5', 6.0', 9.5', 12.0', 14.5', 17.0', 19.5', 22.0' and 25.0'. Pointing 
identifiers are the same as those in Table [T] Extracted point sources are 
shown circled in pink. 




Figure 2. Locations of galaxy clusters in th e galaxy supercluster number 
154 in the catalogue of lEinasto et alj 420011) . shown on RASS hard band 
data. The pointings for Abell 2029 are shown at the bottom of the image. 
Abell 2066 is omitted as it lies 6 degrees to the South East of Abell 2029. 
Abell 2052 is at a different redshift and so not associated with the other 
clusters. 



2 OBSERVATIONS AND DATA REDUCTION 

Five pointings in a cross formation were obtained during obser- 
vations with Suzaku between 2010 January 28 and 30. A further 
pointing from the public archive in the north west direction was 
taken on 2008 January 8 and was also used. An overlay of these 
pointings is shown in Fig.[T] along with the annuli used in the data 
analysis. The annular regions used lie between rings at radii 2.5', 
6.0', 9.5', 12.0', 14.5', 17.0', 19.5', 22.0' and 25.0'. The low expo- 
sure time of the central pointing is compensated by the high surface 
brightness. 

We use HEASoft version 6.11.1, and the calibration database 
(CALDB) released on 201 1 November 10. The raw unfiltered XIS 
data were cleaned using the AEPIPELINE tool which performs the 
standard screening described in Tables 6.1 and 6.2 of the Suzaku 
ABC guidqj, and in addition to the default screening criteria we 
use C0R>6 GeV to restrict the geomagnetic cut off rigidity of 
the particle background , as is typically done jBautz et al.ll2009l : 
iKawaharada et alj uOld) . Calibration regions and bad columns 
were then removed. The defective region in the XISO data (defec- 
tive due to a micrometeorite hit) was excluded. The observations 
were checked to ensure they weren't contaminated by solar wind 
charge exchange (SWCX) as described in section 1331 and the light 
curves of the observations were checked to ensure that no flaring 
had occured. Following the application of the all of this filtering, 
the good exposure times for each pointing are shown in Table [T] 

Ancillary response functions (ARFs) and redistribution matrix 
files (RMFs) are needed to model the response. T hese were created 
using the XISSIMARFGEN jlshisaki et"Zll2007l) and XISRMFGEN 
ftools, respectively. The emission we detect can be divided into two 
broad components; the cluster emission which varies with position 
on the sky, and the background which is assumed constant on the 
sky. Since the spatial variation is different for these two compo- 
nents, their responses will be different. We model the response for 
the background by assuming a uniform background circular region 
of radius 20' [beyond the size of the XIS field of view (FOV)]. For 
the cluster emission we use as input to XISSIMARFGEN a model of 
the surface brightness profile of the cluster. The model we use is 
the beta model from S98 (their equation 5) obtained using ROSAT 
data in the 0.5-2.0 keV band out to 17', though we note that this 
fit is not appropriate in the central region (r<164") due to possible 
cooling flow regions, and that here the surface brightness is higher 
than the model predicts. This model has also been found to be con- 
sistent with A5CA data in both the 0.7-1.5 keV and 3-10 keV bands 
(Mivoshi et al...2005.) . and we show later in section|6]that it is con- 
sistent with our Suzaku da ta. For the central region (r<164") we 
use the beta model found in lLewis et al.l 12003) using Chandra data 
for the surface brightness model, and match the two beta models at 
r=164". 

It has been found previously JGeorge et al.ll2009l) that the re- 
sults are relatively unaffected by changing the surface brightness 
distribution input to XISSIMARFGEN, so we do not expect possible 
deviations from the j3 model in the outskirts to affect the results 
significantly. Indeed, when the analysis was performed using only 
the S98 /3 model as the surface brightness input to XISSIMARF- 
GEN, (which underestimates the surface brightness of the core) the 
results of the spectral fits were essentially unchanged. 

We use only data from the front illuminated (FI) detectors 
XISO and XIS3 owing to the h igher background le vel of the back il- 
luminated (BI) detector XIS 1 JGeorge et al.l2009l) . We use both the 
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Table 1. Observational parameters of the pointings 



Name 


Obs. ID 


Position 


Total Exposure (ks) 


RA 


Dec (J2000) 


1 


804024010 


Centre 


14.2 


227.7409 


5.7498 


2 


804024020 


NE 


43.4 


227.8527 


6.0108 


3 


804024030 


SW 


33.4 


227.6311 


5.4874 


4 


804024040 


SE 


40.4 


228.0053 


5.6496 


5 


804024050 


NW 


36.0 


227.4808 


5.8608 


6 


802060010 


NW (offset) 


289.2 


227.4644 


6.0238 



3x3 editing mode and the 5x5 editing mode. The spectral fits are 
performed by fitting the data for each detector and editing mode si- 
muhaneously. We use a standard ACDM cosmology with Hq — 70 
km s^^ Mpc^^, Qm — 0.3, flA=OJ. All errors unless otherwise 
stated are at the 1 a level. 

All spectral fits were performed in XSPEC 12.7. Ou using the 
extended C-statistic. 



3 BACKGROUND MODELLING AND SUBTRACTION 

Accurate measurements of the low surface brightness outskirts of 
galaxy clusters requires an accurate understanding of the X-ray 
background and its expected level of variation. The X-ray back- 
ground can be separated into two main distinct components. One is 
the cosmic X-ray background (CXB) which consists of unresolved 
point sources of extragalactic origin, mostly believed to be AGN, 
and which dominates the background in the 2-lOkeV band. The 
second is soft foreground emission from within our galaxy which 
is the dominant contributor to the X-ray background at energies be- 
low 2 keV; this is the result of hot gas within our galaxy and so the 
nature of this background component is strongly dependent on the 
direction of observation. We treat each of these background com- 
ponents separately as follows. 



3.1 The CXB 

The CXB has been studied extensively by most major X-ray satel- 
lites, however measurements of its average flux in the 2-10 keV 
band have varied between differen t observatories, as shown in Fig. 
10 of iDe Luca & Molendil J2004h . We take t he total CXB flu x in 
the 2-10 keV band to be the value from Mor etti et all ( |2009|) us- 
ing Swift data, 2.18 ±0.13 x 10"" erg cm"^ s"^ deg"^, which 
is highly consistent with the value obtained using XMM-Newton 
(2.24 ±0.16 x 10'" erg cm~^ s~^ deg~^ in lDe Luca & Molendil 
|2004 ». These value s are higher than those obtained using ASCA 
in iKu shino et al. 1 20021 ) which are used to mo del the CXB in 
iHoshino et al.l ( 1201 Ol) and I Akamatsu et alj 1201 ll) in their studies 
of Suzaku data for the outskirts of Abell 1413 and Abell 2142 re- 
spectively, which is most likely due to the significant stray light 
contamination to which the ASCA observations were subjected. 

Our Suzaku observations allow us to exclude point sources 
down to a threshold flux of ^exci = 10"^"^ erg cm~^ s"'^ in the 
2-10 keV band (and the point sources excluded are shown in Fig. 
[Hi. To calculate the unresolved CXB level following the removal 
of these point sources, we use the source flux distribution found 
in lMoretti et al . (2003J (obtained by fitting to a catalogue of point 



sources extending down to a flux of 2 x 10"^" erg cm~^ s~^) 
to calculate the flux contribution of the resolved sources. We do 
this instead of integrating the source flux distribution from to the 



threshold flux because the source flux distribut i on is u ncertain be- 
low the minimum flux studied in lMoretti et alj (12003'). and indeed 
when integrated over all ener gies leads to an und erestimate of the 
total CXB level (see Fig. 5 of lMoretti et alj2003h . 

The unresolved CXB level (in ergs cm~^ s~^ deg'^) is there- 
fore given by: 



FcxB = 2.18 ±0.13 X 10" 



(§)=<"- 



(1) 



which, using the two powerlaw best fit for N(S) given 
in iMoretti et alj ( |2003|) gives an unresolved 2-10 keV flux of 
1.87±0.13 X 10"" ergcm"^ s'^ deg'^ 

The expected fluctuations in the CXB for a given solid angle, 
SI, due to variations in the number of unresolved point sources be- 
low th e threshold flux in the area of observation can be calculated 
using ( JBautz et alj2009l) : 



O-CXB = (1/0) 



"'©-'- 



(2) 



where we see that the variance decreases with increasing solid 
angle observed (fJ) and with decreasing threshold flux (Sexci)- 

The expected fluctuations in the CXB for each extraction re- 
gion examined are shown in table |2] These provide the ranges 
which must be investigated to understand the effect of the uncer- 
tainty of the CXB level on our measurements in the outskirts. These 
ranges are consistent with those obtained by scaling the fluctuations 
observed with Ginga using the method of Hoshino et al. (2010Q, 
and also with the observed variance described in iSimionescu et al.l 
(20ii) using the Suzaku background pointings for the Perseus clus- 
ter. 



3.2 The soft galactic foreground 

Abell 2029 is located well out of the galactic plane, meaning that 
the absorbing hydrogen colu mn densi ty is low [the mean Leiden 
Argentine Bonn (LAB) survey dKalberla et alj2005h value is 3.26 x 
10^'' cm~^] and that possible variations of the column density 
across the cluster have a negligible affect on the soft band emis- 
sion. 

Abell 2029 is however located near the North Polar Spur 
(NPS), which is a region of complex soft band emission attributed 
to part of a hot interstellar bubble formed from supernova explo- 
sions and the outflows of young, hot stars. Previous studies of Abell 
2029 and its surrounding region have required extra background 
components to model the NPS emission, in addition to the back- 
ground level associated wit h normal galactic em ission. Studies of 
Abell 20 29 with Chandra (fmhlinin et al.ll200^ and with XMM- 
Newton JBourdin & Mazzottal 120081) have used a soft foreground 
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Table 2. CXB fluctuations in each annulus for each pointing in the 2-10 keV band (10 ^^ erg s ^ cm ^ deg ^) 



0.0'-2.5' 2.5'-6.0' 6.0'-9.5' 9.5'-12.0' 12.0'-14.5' 14.5'-17.0' 17.0'-19.5' 19.5'-22.0' 22.0'-25.0' 



Position 



Centre 

NE 

SW 

SE 

NW 

NW (offset) 



6.8 



3.1 



2.7 



6.1 


5.0 


4.6 


4.5 


4.5 


4.6 


5.6 


5.0 


4.6 


4.5 


4.6 


4.9 


5.5 


4.6 


5.0 


4.7 


5.2 


6.4 


6.6 


6.5 


4.9 


5.1 


5.3 


5.0 


- 


- 


8.2 


5.6 


5.7 


5.1 



model consisting of two unabsorbed MEKAL components, one at 
0.46 ± 0.07 keV and the other at 0.22 ± 0.05 keV. 

We use the region outside 25' from the north-western offset 
Suzaku pointing to fit for the soft thermal background components, 
and find that that the best fit (shown in Appendix A) is achieved 
using a model consisting of two unabsorbed APEC components 



with temperatures at 0.53l^:^g keV and Q.IQZq^I keV, each of 
which has its metallicity fix ed to 1 Z^. This is c onsistent with 
the foreground model used in lVikhlinin et al.l ( 120050 described ear- 
lier, which themselves were found to be consistent with previous 
Chandra measurements of the NFS region. To investigate the spa- 
tial variation of these background components, we use the ROSAT 
PSPC pointing of Abell 2029 (observation identifier: rp800249) 
and extract four background regions (in an annulus between 30' to 
50') to the northwest, sorthwest, northeast and southeast, and fit the 
spectra in the 0.4-2. OkeV bands after subtracting the particle back- 
ground as modelled by the ftool PCPARPHA. ROSAT ARFs were 
created using the ftool PCARF and standard filtering procedures 
were used as described in the document 'ROSAT data analysis us- 
ing xselect and f tools' Q We find that the spectra are completely 
consistent with the two APEC component model found using the 
Suzaku region to the north-west, and table|3]shows the spatial vari- 
ations in the normalisations of the two background components, 
which we use as a indicator of the extent of expected spatial varia- 
tions of these soft background components. The fits to the ROSAT 
data are also shown in Appendix A. These variations will be used 
later, along with the variations expected in the CXB from table|2] to 
determine the affect on our fits of our uncertainty in the background 
model. 

When fitting for the soft foreground the CXB is modelled as an 
absorbed powerlaw of index 1.4 and its normalisation is allowed to 
be a free parameter In all cases the CXB normalisation is consistent 
with the expected level given the level of point source removal. 



3.3 Non X-ray Background 

We model the non X-ray background (NXB) resulting from inter- 
actions with charged particles using the XISNXBGEN ftool, using a 
database of measurements taken in an interval spanning 150 days 
before and 150 days after the observations, and subtract this from 
the spectra. 



3.4 Stray light modelling. 

Due to broad PSF each annulus contains contributions from 
neighbouring annuli. In addition, the precollimator is not com- 
pletely successful in preventing light following undesired paths 



heasarc.gsfc.nasa.gov/docs/rosat/ros_xselect_guide/xselect_ftools.html 



(such as back reflection and secondary reflection, as described in 
ISerlemitsos et al.l 120071) , and so stray light enters each pointing 
from the bright core. 

To understand the magnitude of this effect we use the image 
of each annulus (taken to be represented by the beta model of S98 
which we show later fits our Suzaku data well) as input to the ray 
tracing simulator XISSIM and simulate a 2x10^ photon exposure 
for each annulus, and then measure the fraction of photons from 
each annulus which are spread into the other annuli. The resulting 
distribution of photons is shown in table|4] We see that the majority 
of the flux detected in each annulus originated from that annulus, 
and that in the outskirts the majority of the contamination originates 
from the nearest innermost annulus which has similar properties to 
the annulus in question. 

We model emission from the central 6.0' (from which the ma- 
jority of the cluster flux originates) which spreads into the arms of 
the cross using XISSIM, using the spectra and image of the core 
as inputs. We model this emission and include it as a background 
component, and investigate the affect of accounting for this stray 
light (i.e. the systematic error caused by the stray light) on our fits 
in section|43] 



3.5 Solar Wind Charge Exchange 

To ensure our observations are not significantly contaminated due 
to solar wind charge exchange (SWCX) emission we use data 
from the Solar Wind Ion Composition Spectrometer (SWICS) on 
the Advanced C omposition Explorer (ACE) spacecraft. We follow 
ISnowden et alj 112004). which assumes that the SWCX is negligible 
when the ratio of O^^ to O"''' is below ~ 0.2 (this criterion is also 
used in Humphrev et al. 2011), and we show in figure |3]that only 
very small time intervals need to be excluded because of this, and 
that O^^ 10^^ does not increase substantially above 0.2 during the 
times of observation. 



4 MODELLING CLUSTER EMISSION 

We fix the X-ray background model using the CXB level deter- 
mined in section [3TT] following point source subtraction and model 
the CXB as an absorbed powerlaw of index 1 .4, added to the galac- 
tic components (two APEC components) found in section lT2l bv fit- 
ting to the outer annuli of the north western offset pointing (which 
are found to be consistent with ROSAT results for a background 
annulus between 30'-50' around the cluster). Later we will inves- 
tigate the systematic error of fixing the background by varying its 
parameters within the confidence regions we determined in section 



m 
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Table 3. Soft foreground fluctuations from ROSAT PSPC and Suzaku data. The units of the APEC normalisations are 10-^4(4^)-i£)-^(X _l_ 
z)— ^ f nenndV, where Da is the angular size distance (cm), and Ue and nif are the electron and hydrogen densities (cm-3) respectively, and these 
values are scaled for a circular area of sky of 20' radius (1257 arcmin'^). 



Position 




0.53 keV APEC norm 


0.20 keV APEC norm 


ROSAT NE 

ROSAT SE 

ROSAT NW 

ROSAT SW 

SUZAKU NW offset (25' 

SUZAKU NW offset (29' 


-29') 
-34') 


5.5+°-^ X 10-4 

e.itg;? X 10-4 

6.1+2-^ X 10-4 
6.9+0-^ X 10-4 
6.8+|;-;^ X 10-4 
6.7+0-5 ^ 10-4 


3.0+0^ X 10-3 

3.61UI X 10-3 

2.9+01 X 10-3 
4.3+0^3 ^ 10-3 

3.5t^j X 10-3 
3.4+°;3 X 10-3 



Table 4. Percentage contribution of flux in the rows' annulus from the columns' annulus due to PSF spreading and stray light. 





0.0'-2.5' 


2.5'-6.0' 


6.0'-9.5' 


9.5'-12.0' 


12.0'-14.5' 


14.5'-17.0' 


17.0'-19.5' 


19.5'-22.0' 


0.0'-2.5' 


95.9 


3.9 


0.05 


0.002 


0.0002 


0.0004 


0.0004 


0.0001 


2.5'-6.0' 


36.2 


60.2 


3.5 


0.04 


0.009 


0.004 


0.003 


0.002 


6.0'-9.5' 


14.5 


19.8 


64.4 


1.1 


0.07 


0.02 


0.02 


0.01 


9.5'-12.0' 


6.7 


3.5 


15.9 


63.7 


9.1 


0.7 


0.2 


0.05 


12.0'-14.5' 


7.6 


3.5 


4.0 


13.8 


62.5 


7.3 


0.8 


0.2 


14.5'-17.0' 


4.3 


1.6 


2.3 


1.9 


15.4 


62.8 


11.1 


0.6 


17.0'-19.5' 


1.3 


1.0 


1.3 


1.1 


2.0 


19.8 


65.8 


7.4 


19.5'-22.0' 


5.2 


1.2 


1.0 


0.7 


1.1 


3.0 


19.8 


65.1 



4.1 Projected fits 

First of all we do projected fits, modelling each annulus as an APEC 
component absorbed by a column density of 3.26x10^" cm^^ 
(from the LAB survey) using the PHABS absorption model. We 
fit to 0.7-7.0 keV band and for the central 3 annuli we allow the 
metallicity to be a free parameter. In the 0.0'-2.5', 2.5'-6.0' and 
6.0'-9.5' annuli the best fit projected metallicities are 0.45lo o^ ^e 



1.03 ■^0' 



0-3lQQg Zq and 0.2_oi Zq, which are r easonably consistent with 
the metallicities found in lVikhlinin et al.l |J2005). In the annuli out- 
side 9.5' the metallicity is fixed t o 0. 1 Zq as it is not c onstrained 
in fits (this is the value obtained in lVikhlinin et al.ll2005l in the out- 
skirts). We find that varying the metallicity in the outskirts through 
a reasonable range (from Zq to 0.3 Zq) has a negligible affect 
on the fitting parameters. We fix the redshift to 0.0767. We tie the 
temperatures of the outer 3 annuli due to the low number of counts 
in the outermost annuli, but let their normalisations be free between 
them. 

We initially examined each direction separately and searched 
for significant azimuthal variations which would make averaging 
over certain directions inappropriate. We find two main sources 
of asymmetry: firstly the temperature is lower in the south east- 
ern direction than in the other directions which are consistent with 
one another. This is shown in Fig.|6](the deprojected temperatures 
shown there are consistent within errors w ith the projected t em- 
peratures, which has also been observed in iBautz et al .1120091 and 
lAkamatsu et al.l201 1 ). This lower temperature to the south east co- 
incides with a depression in temperature se en in the temperature 
map produced using XMM-Newton data in iBourdin & Mazzottal 
uOOa to the south east, suggesting a disturbance of the ICM in this 
direction extending from the core of the cluster all the way to the 
outskirts, possibly due to the accretion of small galaxy groups in 
this direction. 

Secondly the projected density is highest to the north and flat- 



tens off (see Fig. |4]l and this is the only direction in which we de- 
tect statistically significant emission beyond 22'. This indicates the 
presence of excess emission above the azimuthal average which 
may indicate a filament structure to the north connecting Abell 
2029 and Abell 2033, or emission from Abell 2033 entering the 
field of view. We investigate this in depth in section|6] 



4.2 Deprojected fits 

We perform two deprojection runs. The first averages over all of the 
directions except for the northern direction (where excess emission 
exists, rendering deprojection impossible in this direction), and the 
south eastern direction (which has a lower projected temperature 
and entropy than the other directions). In the second the south east- 
ern direction is deprojected individually due to its lower tempera- 
ture and entropy. 

To achieve the deprojection we fit all nine annuli simulta- 
neously, with each annulus modelled as the superposition of the 
emission from the shell associated with the annulus and the emis- 
sion projected into that annulus from shells exterior to it under 
the assumption of spherical symmetry (thus emulating the PRO- 
JC T mixing model by add ing multiple APEC components as is done 
in lHumphrev et al.ll201 iL which allows data from multiple Suzaku 
instruments and pointings to be fitted simultaneously). The tem- 
perature is tied for the 3 annuli between 17' and 25' and for the 
3 annuli between 9.5' and 17' due to the limited spectral quality 
of the data. The resulting deprojected temperatures are consistent 
with the projected ones. The deprojected temperature, density and 
entropy (S=kT/ne ) profiles are shown in Fig.|6]for the azimuthal 
average and the south eastern direction separately. 

The lower temperature in the south eastern direction causes 
the entropy profile to flatten indicating significant deviation 
from hydrostatic equilibrium. The temperature profile to the 
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Figure 3. 0^+/0®+ ratio as measured by the ACE spacecraft. SWCX is assumed to be negligible when 0^"'"/0®+ is below ^ 0.2, and regions above this 
level are excluded. Left is for the observation taken in 2008 (802060010), and right is for the time region of the observations in 2010 (804024010, 804024020, 
804024030, 804024040, 804024050). The day of observation time axis label is the day of the year on which the observation took place. 



south east maps on well to the temperature profile found in 
iBourdin & Mazzottal 12008) for the low temperature region they 
find to the south east (within 8 arcmins of the core), and is overplot- 
ted in blue in Fig.|6] The lower entropy profile from XMM-Newton 
in the central 8 arcmins also maps on well to the flatter and lower 
entropy we find in the outskirts to the south east, suggesting that the 
ICM has been disturbed from equilibrium from the outskirts to near 
the core, possibly due to the accret ion of small galaxy groups a long 
this direction (as was suggested in lBourdin & Mazzottall2008fi . 

Due to Suzaku 's large PSF it is not possible to divide the in- 
ner two annuli into sectors corresponding to the arms of the cross, 
but when the 3rd annulus (between 6.0'-9.5') is divided we find a 
statistically significant drop in the temperature below the azimuthal 
average, which is consistent with the XMM-Newton temperature in 
the same region, as shown in the left column of Fig.|6] 

For the azimuthally averaged fits the entropy profile rises with 
increasing radius , but less steeply than the r^'^ relation derived in 
IVoit et al.l ( 120050 assuming purely gravitational hierarchical struc- 
ture formation (and neglecting non gravitational processes such as 
feedback, cooling and star formation). 

The temperature profile for all directions neglecting the south 
east is consistent with what has previously been found using Chan- 
dra, Swift and XMM-Newton, as shown in Fig.|5]and Fig.|6] The de- 
projected density profile in all directions is consistent with the ex- 
trapolated best fit model to Chandra data in (Vikhlinin et al. 2006) 
when systematic errors in the background modelling are taken into 
account, as shown in Fig.|6] 
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Figure 4. Comparing projected densities. Grey is the azimuthal average 
excluding the north. Blue is for the low temperature south eastern direction. 
Red is for the northern direction, showing a significant excess. Solid lines 
show 1 a systematic errors calculated as desciibed in section 1431 



4.3 Systematic errors 

We need to understand the contribution of the uncertainty of the 
background model to the error budget. The normalisations of the 
CXB and the galactic components, (and the temperatures of the 
galactic components) are only known to an accuracy described in 



section [3] so we need to see the effect on our fits by exploring 
the possible background parameter space. Typically (for example 
George et al. 2009 and Akamatsu et al. 20 11), these systematic er- 
rors have been quantified by varying only one parameter at a time 
and leaving the others unchanged, which underestimates the actual 
uncertainty. 
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Figure 5. Comparing Suzaku temperature profiles with projected tempera- 
tu res obtained with Swift in|Moretti et al. (2011) (green) and Chandra (red) 
in Ivikhlinin et al.1 (2009). Grey crosses are azimuthally averaged Suzaku 
temperatures, (exluding the north and the south east). The black dashed line 
is the deprojected temperature from Chandra in lvikhlinin et alj (200g) and 
the sohd black line shows the best fit profile from 6' to 22' to the depro- 
jected Suzaku temperatures. 



Here we compute the systematic uncertainty of the back- 
ground model on our deprojected fits by producing 10000 random 
realisations of the background parameters [the CXB normalisation, 
the 0.53 and 0.20 keV galactic APEC component normalisations 
and temperatures, a nd the NXB leve l (which has an uncertainty of 
± 3.6 percent from lTawa et al.ll2008i) distributed by their expected 
variance and performing the deprojected fit for each. In this way 
the systematic uncertainty in the background is folded into the de- 
projection to give a complete propagation of the uncertainty in the 
background to the uncertainty in the spectral fits. The one sigma 
systematic errors on the temperature, density and entropy profiles 
resulting from this process are shown for the deprojected fits in Fig. 
|6] For the temperature the systematic error is roughly of the same 
order as the statistical error (except in the outermost bin), but for 
the density it is the dominant source of uncertainty in the outer 4 
annuli. 

The systematic errors are also calculated for the projected den- 
sity and shown in Fig.|4l showing that the excess to the north is still 
statisically significant. Note that because these are projected fits the 
errors are smaller than for the deprojected fits (because for the de- 
projected fits the error on the outer annuli is propagated to the error 
on the inner annuli). 

We then perform all the deprojected fits with the stray light 
level from the central 6' (as simulated by XISSIM) subtracted, and 
plot as the green line in Fig.|6]the best fit results with this stray light 
subtracted. In all cases the effect of taking into account the stray 
light is within the systematic error of the background modelling. 

We vary the contamination level on the optical blocking filter 
(OBF) by ±10 percent as described in lAkamatsu et al.l ( 1201 ih . us- 
ing the Suzaku ftool XISCONTAMICALC to modifty the ARFS, and 
find this to have a negligible affect on the fits. 



5 MASS ANALYSIS 

Assuming hydrostatic equilibrium an d spherical syrnmetry the total 
mass interior to radius r is given by (.Vikhlinin et al.ll2006n 



M(< r) = -3.68 X lO"M0r(r-)r 



d In nn dlnT^, 

+ — 1 (3) 



din r 



din r 



where T is in keV and r is in Mpc. We find best fit tem- 
perature and density profiles by using Monte Carlo methods with 
lOOOO trials, and use these profiles in equation|3]to yield the grav- 
itational mass profile and its 1 a error. The density profile is fit 
as a beta model of the form nH=no(l-l-(r/rcn)^)~ , and the tem- 
perature profile is fit with the function T=To(l+r/rcT)" as used in 
lAkamatsu et al.l ( I2OIII) and lBums et ahl ( |20I(]|) . We do this for the 
azimuthally averaged profiles shown in the left hand column of Fig. 
[6](thus excluding the north due to its excess emission and the south 
east due to its lower temperature and entropy which means that it 
is likely out of hydrostatic equilibrium). From this we then calcu- 
late an estimate of T200 (the radius within which the mean density 
is 200 times the critical density of the universe), which is widely 
used as an approximation to the virial radius. The cumulative gas 
mass is calculated using the number density of the ICM derived 
from the APEC normalisation and assuming spherical symmetry, 
allowing the gas mass fraction, fgas, as a function of radius to be 
derived. 

The slig htly shallower and lowe r temperature profile we find 
compared to Ivikhlinin et al.l (|2006), together with the slightly 
higher density in the outskirts causes our total mass estimate to 
be lower, and gas mass estimate to be higher, making our mea- 
surement of the gas mass fraction in the outskirts higher than those 
presented in Vikhlinin et al. (2006). This consequently means that 
our estimate of r2oo is smaller. We find M5oo=7.2lVj, x 10^* 
M0, which is consistent with the IVikhlinin et alj (200g) value of 
M5oo=8.01±0.74 X 10" Mq. We find r2oo=1.92lg:iJ Mpc (cor- 
responding to 22.0l:i;4 arcmins), with M2oo=8.0li;5 x 10^* M©, 
which is smaller than the value obtained using the scaling relation 
of lAmaudetal.l ( l2005h : 



M200/MQ = 5.74 ± 0.3 X 10" [^j-^J /E{z)^% 

E{z) = [^^{l + zf + l-n^\ (5) 

(6) 

which yields M200 =12.lto;g x 10 " Mq andr2oo=25.3to.6 
using the mean spectroscopic te mperature o f 8.47±0.09 keV from 
Chandra obtained in .Vikhlinin et alj J200Q) . Part of this discrep- 
ancy is, as described in lBautz et alj ( 120091) (which similarly finds a 
lower mass and r2oo than derived using Chandra results and scaling 
relations), likely due to the fact that Chandra reports temperatures 
for hot clusters which are systematically higher than those from 
XMM-Newton. This can be seen by comparing the XMM-Newton 
temperature points in Fig. |6] with the Chandra points in Fig. [5] 
and it is believed to be due to the uncertainty in the contamination 
on the Chandra mirrors. Our Suzaku results suggest a lower emis- 
sion averaged spectroscopic temperature of 7.5±0.1 keV (which is 
consistent with the average temp erature of 7.45±0. 19 found with 
ASC A in iFukazawa et alj |2004| and is more consistent with the 
XMM-Newton temperatures), which when used in equation|4] gives 
M2oo=10.i1o;6 X 10^" Mq andr2oo=23.8lo;5 arcmins, which are 
consistent with the values derived from our Suzaku fits. 

The gas mass fraction is shown in Fig. |7] which rises with ra- 
dius and reaches the cosmic mean baryon fraction (0.167, obtained 
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Figure 6. Deprojected temperature, density and entropy profiles. Left column is for the azimuthal average (excluding the north which has excess emission 
and the south east which has a lower temperature and density). Right column is for the south eastem direction which has a lower outskirts temperature and 
flatter entropy profile. Solid black lines show one sigma systematic errors from uncertainties in background modelling. The solid green line shows the e f fect of 
subtracting the stray light emission from the core. The red points overplotted on the right are the XMM-Newton temperatures from Bourdin & Mazzottaj fcOOSb 
found for the region away from the dep ression in temperature to t he south east. The blue points for the temperature and entropy on the left are the for the cold 
depression to the south east found in lBourdin & Mazzottal (2008) and correlate well with the Suzaku values in the outskirts to the south east. For the density 
plots the solid red line i s the best fi t deproiected density profile found in Vikhlinin et al. (2006) using Chandra data. For the entropy plots the solid black line 
is the r^'^ relation from lVoit et al.l ((2005) scaled to pass through the central three points of our Suzaku data on the left. The same profile is shown to the right 
to emphasise the lower entropy in the south east. 



from CMB data in lKomatsu et alj|201 ll ) near our calculated value 
of the virial radius. 



6 EXCESS TO THE NORTH. 

The projected density profile shows an excess to the North, sug- 
gesting the existence of a filament connecting Abell 2029 to Abell 
2033, or contaminating emission from Abell 2033 entering the field 
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Figure 7. Gas mass fraction as a function of radius for Abell 2029. Black 
line shows the best fit and the green lines show the 1 cr error ranges. The 
red h orizontal line shows the cosmic mean baryon fraction as determined in 
iKoma tsu et al. 1 201 1). The solid blue vertical line shows our measurement 
of r2oo and the vertical blue dotted lines show the one sigma uncertainty in 
r200- 



of view. The existe nce of this stmcture has also been claimed in 
lEckeit et alj ( 120 llh who studied azimuthal variations in the sur- 
face brightness profiles of 3 1 galaxy clusters using archival ROS AT 
PSPC data, and found a significant excess connecting Abell 2029 
and Abell 2033. 

To check the presence of this excess we used a ROSAT PSPC 
pointing of the Abell 2029 and Abell 2033 system (rp800249) and 
examined the surface brightness profiles in the directions corre- 
sponding to the arms of the cross (using the same regions and annuli 
binning as used for our Siizaku pointings). The ROSAT PSPC data 
were reduced using the ROSAT E xtended Source Analysis Soft- 
w are (ESAS, ISnowden et al.|[l994l) . We follow the same method 
as lEckert et al.l ( 120111) to obtain a background subtracted, exposure 
corrected counts image for the R47 band (corresponding to 0.44- 
2.04keV). The resulting image is shown in Fig. [S] and there is a 
clear bright region connecting Abell 2029 to Abell 2033. 

The resulting background subtracted projected ROSAT sur- 
face brightness profile is shown in Fig. |9] where the error in the 
background has been added in quadrature to the error bars. We 
clearly see an excess to the north (red) deviating from the beta law 
decrease of the azimuthal average of the other directions. This is 
the only statistically significant deviation from azimuthal symme- 
try in the surface brightness profiles. We also show the correspond- 
ing vignetting corrected and background subtracted projected sur- 
face brightness profile obtained with Suzaku in the 0.7-2. OkeV band 
in Fig. |9] (where vignetting correction and background subtr action 
were performed following the method of ISautz et al.ll2009h . also 
showing the excess in the northern direction. The overplotted beta 
model is that obtained in S98 (which used different data reduction 
methods) as the best fit to the azimuthally averaged profile between 
164 arcsecs and 1000 arcsecs. 

However another possibility is that the excess is caused due to 




Figure 8. Exposure corrected and background subtracted ROSAT PSPC 
image of the Abell 2029/2033 system smoothed with a Gaussian kernel, 
showing indications of excess emission to the north of Abell 2029 between 
the two clusters. The virial radii of the two clusters are shown by the blue 
circles. 



overlapping emission from the outskirts of Abell 2029 and Abell 
2033, with no need for excess emission from a filament. Using the 
scaling relation, equation [4] for r2oo we find r2oo for Abell 2033 
to be 1.5 Mpc (17.0' ) (using the mean tem perature of 3.84 keV 
from ASCA found in iFukazawa et al.ll2004l) . We found in section 
[5]that for Abell 2029 our Suzaku results imply r2oo=22.0!ti'4 ar- 
cinins. These virial radii are overplotted as the blue circles on Fig. 
[8] and we see that the excess emission correlates well with the re- 
gion where we would expect to detect emission from the ICM of 
the outskirts of Abell 2033. 



6.1 Spectral analysis of northern excess 

The difference in redshift between Abell 2029 (z=0.0767) and 
Abell 2033 (z=0.081) corresponds to a distance of ~ 20 Mpc (ne- 
glecting the effect of possible peculiar velocities), meaning that if 
a filament exists between the two then we are seeing it along its 
length, increasing its surface brightness due to this preferential ge- 
ometry. Such a preferential geometry has already been exploited 
to study a potentia l filament between Abell 222 and Abell 223 
dWemer et al.ll2008h which is viewed along its length, however as 
that system is located at a much higher redshift (z=0.21) than the 
Abell 2029 and Abell 2033 system its surface brightness and area 
of extent are both lower. 

To fit for the projected properties of the northern excess, we 
perform a fit to the 19.5'-25.0' region to the north, and simultane- 
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Figure 9. Left: Exposure corrected and background subtracted ROSAT PSPC surface brightness profile for the same regions analysed with Siizaku of the Abell 
2029/2033 system in the 0.4-2.0 keV band. RightiVignetting corrected and background subtracted Suzaku surface brightness profile in the 0.7-2.0 keV band, 
averaged over the two front illuminated detectors (XISO and XIS3). In both panels the black points are an azimuthal average over all directions except the 
north and the red points are for the northern direction only, showing a significant excess to the North with both satellites. The black solid line is the beta model 
found in S98 for the azimuthally averaged surface brightness profile between 164 arcsecs and 1000 arcsecs (16.7 arcmin). 



ously fit the same regions in the other 4 directions which are tied 
to get an azimuthal average, and use this as a background model. 
In this way we fit the excess over the azimuthal average level to 
the north in the 19.5'-25.0' region as an absorbed APEC compo- 
nent. The excess of the spectrum over the azimuthal average of 
the other directions is shown in Fig. [To] The normalisation of the 
excess APEC component, taking into account systematic errors in 
the background determination is l.Sto's ^ 10 ~'^ cm~*, and the 
temperature is 3.6l]^'2 keV). This temperature appears too high for 
the excess to be caused by a WHIM filament connecting the two 
clusters (which would be expected to have a temperature around 1 
keV as t his gas wou ld not ha ve been shock heated due to accre- 
tion yet j Werner et a l. 2008, D ave et akluOOlh . The temperature is 
more consistent with this emission originating from the outskirts of 
Abell 2033. Assuming that the excess emission completely origi- 
nates from Abell 2033 and assuming spherical symmetry, the emis- 
sion norm of the excess corresponds to a density in the range 1 .0- 
2.0 X 10 ""' cm~^ in the outskirts of Abell 2033, which seems 
reasonable. 

Figures|4]and|9]show that the excess to the northern direction 
is present from a radius of 12' outwards, but it only seems rea- 
sonable to attribute the excess at radii beyond 19.5' as being due 
to overlapping emission from Abell 2033. The excess between 12' 
and 19.5' could have a number of causes. One is that the density 
to the north could genuinely be higher than that in other directions. 
Another is that we may be detecting emission from gas in a fila- 
mentary structure connecting Abell 2029 and Abell 2033. 



6.2 Comparing to the spatial galaxy distribution 

To investigate the nature of the excess emission between Abell 
2029 and Abell 2033 further, we constructed the density field of 




Energy (keV) 

Figure 10. Spectrum of excess emission to the north. Points show spectrum 
of northern excess region. Line shows the azimuthal average spectrum fit 
for all of the other directions. 



galaxies around Abell 2029 using data from the Sloan Digital Sky 
Survey (SDSS). Taking the photometric catalogue from SDSS Data 
Release 7 (Abazajian et al. 2009), we select galaxies with mag- 
nitude i < 2L0 and a photometric redshift within a range of 
dz = 0.032 centred on the redshift of Abell 2029, Zc = 0.0767. 
This range corresponds to the median photometric redshift uncer- 
tainty for sources of this magnitude at 2 < O.L The general fea- 
tures of the density field are not highly sensitive to the redshift 
range used, and the spectroscopic catalogue produces similar re- 
sults though with a lower sampling density. Though the range ex- 
ceeds the typical velocity dispersion of a cluster, we can easily 
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Figure 11. This shows the distribution of galaxies with magnitude i<21.0 
and photometric redshift \z — Zc\ <0.016 using SDSS data. We define 
the overdensity as (p — {p))/(p) where (p) is the mean galaxy density 
computed on the 16 deg^ region centred on Abell 2029, and smooth the 
overdensity plot using a Gaussian kernel with a FWHM of 500 kpc at the 
redshift of Abell 2029. 



identify peaks in the density field associated with known clusters 
nearby, as shown in Fig. \TT\ We plot the overdensity of galaxies 
smoothed by a Gaussian with FWHM= 500 kpc. 

There is a clear excess in the galaxy density to the north which 
appears to originate from Abell 2033, consistent with our hypoth- 
esis that we are seeing the overlap of the virial radii of Abell 2029 
and Abell 2033. 

The foreground object NSC J15 1109+052006 is very faint in 
X-rays (it is visible as an extended source to the south in Fig. [8} 
and lies out of the field of view of the south west pointing. We 
simulated the expected contribution from this source in the south 
west pointing using XISSIM and found it to be negligible. 



7 SUMMARY 

We explore the ICM of the galaxy cluster Abell 2029 to higher ra- 
dius than before, with no statistically significant emission detected 
beyond 22' (except for the northern pointing between Abell 2029 
and Abell 2033), which is 1 .9Mpc, and is equal to our measurement 
of r2oo • Excess emission is found to the north between Abell 2029 
and Abell 2033 which is detected in both our Siizaku observations 
and in ROSAT PSPC observations, and the most consistent expla- 
nation for this excess emission is that it originates from Abell 2033 
and we are seeing the overlap of the outskirts of the two clusters to 
the north. 

We detect a cold feature to the SE extending out to the edge of 
the detected cluster (22') where the ICM is significantly colder than 
in the other directions, consi stent with the XMM-Newton findings of 
iBourdin & Mazzottal ( 120081) which found a temperature depression 
to the south east within the central 8'. The substantial point spread 
function does not allow this cold feature to be traced into the central 
2 annuli of the Suzaku pointing. The lower temperature causes the 
entropy profile to be lower and to flatten, indicating that the ICM 
in this direction is significantly out of hydrostatic equilibrium. This 



may be the result of the accretion of galaxy groups along the SE 
direction, which has disturbed the ICM. 

Away from this cold feature, the azimuthally averaged tem- 
peratures and densities are consistent with previous findings, with 
the entropy profile rising in the outskirts but not as steeply as 
the r^'^ powerlaw predicted due to pure gravitational hierarchi- 
cal structure formation, a result which has also been found for 
other clusters [Abell 17 95, JBautz et al1l2009l). the Perseus clus ter 
JSimionescu et al.l20Ill) . and Abell 2142 JAkamatsu et al.ll20lll) 1. 
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Figure Al. Fitting to the Suzaku background regions in tlie NW offset pointing. Left sliows tlie fit to tlie 25'-29' annulus, and riglit sliows tlie fit to tlie 29' 
to 34' annulus using. Overplotted in green in the NXB level. Note that the areas of extraction for the two spectra are different, thus explaining the different 
normalised count rates on the y axis. 
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Figure A2. Fitting to ROSAT PSPC background regions in an annulus between 30' and 50' using the same model found using the Suzaku background pointings. 
Top left is for the NW, top right is for the NE. bottom left is for the SE and bottom right is for the SW. These fits give the ranges in the normalisations of 
the background parameters shown in table |3]which are used to understand the expected spatial variation of the soft X-ray background. Note that the areas of 
extraction for the four spectra are different, thus explaining the different normalised count rates on the y axis. 
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APPENDIX B: SPECTRAL FITTING 
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Figure Bl. Spectral fitting to the XISO and XIS3 data simultaneously. The spectra are for the azimuthally averaged fit (excluding the north and the south 
east) used to obtain the results in the left column of Fig. |6] but here we show the projected fits for simplicity. In reading order the spectra correspond to 
increasing radii annuh, which are between 0.0'-2.5', 2.5'-6.0', 6.0'-9.5', 9.5'-12.0', 12.0'-14.5', 14.5'-17.0', 17.0'-19.5' and 19.5'-22.0'. The central 3 annuh 
are integrated across all azimuth and are from the central pointing only. The red lines through the points represent the best fits (background plus cluster 
emission), while the lower red line shows the excess over the background level (shown as the black solid line) modelled in these projected fits as an absorbed 
APEC component. The NXB level is overplotted in green. 
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